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Abstract 

Assuming that the nucleus can be treated as a perfect fluid we study the conditions for the 
formation and propagation of Korteweg-de Vries (KdV) solitons in nuclear matter. The KdV 
equation is obtained from the Euler and continuity equations in nonrelativistic hydrodynamics. 
The existence of these solitons depends on the nuclear equation of state, which, in our approach, 
comes from well known relativistic mean field models. We reexamine early works on nuclear 
solitons, replacing the old equations of state by new ones, based on QHD and on its variants. Our 
analysis suggests that KdV solitons may indeed be formed in the nucleus with a width which, in 
some cases, can be smaller than one fermi. 
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I. INTRODUCTION 


Understanding the nucleus and its properties in terms of QCD is a very ambitious goal 
and we are still far from it. Significant progress has been achieved in this direction with 


the use of chiral perturbation theory Q|. On the other hand, relativistic mean held (RMF) 
models have been sucessfully'used to describe the bulk properties of nuclear matter, hnite 
„uc,ei and denae .a. A.don,. ... a. on,. edeCve .node,a. tde. ane a,.np,e, 

depend only on few parameters and are the ideal tool to test certain aspects of nuclear 
matter. 

In this work we shall use the relativistic mean held model discussed in P| to investigate 
the formation and propagation of Korteveg - de Vries solitons in nuclear matter. Assuming 
that nuclear matter behaves like an ideal huid we wish to know whether small baryon density 
perturbations can propagate without dissipation through the medium. If this turns out to be 
case, these pulses might be responsible for interesting phenomena. One of them would be the 
apparent nuclear transparency in proton nucleus collisions at low energies. The impinging 
proton would be absorbed by the huid and turned into a density pulse. It would then be able 
to Uavers^the target nucleus emerging on the other side. This idea has been hrst advanced 
the experimental evidence of this “nuclear transparency” was discussed. In 


LU LiClVCiDC LJ 

in 1^. In 


the formalism was improved with the inclusion of viscosity and three-dimensional expansion. 


In j9| the ehect of the nuclear surface was considered. One important lesson from all these 
works is tlmt^l^ formation of solitons depends in a crucial way on the equation of state 
(EOS). In 


the equation of state employed was borrowed from Q| and it was very 
phenomenological, containing no information about the microscopic nucleon dynamics. 

The experimental signature of the nuclear soliton is very difficult to measure, since the 
outcoming nucleon emerges in a direction very close to the beam. Probably because of this 
reason the subject was left aside for almost ten years. Soliton formation in nuclear matter 


developed 


in Q 


Q. 


was then revived in [11|, where a different equation of state was considered. This EOS was 


and was based on the nucleon-nucleon microscopic interaction given by the 


Skyrme force. The main Ending of ^ was that the nuclear soliton is not a compression 
pulse, as thought before. Instead, it is a rarefaction pulse due to a sign inversion appearing 
because of the use of Skyrme interactions. 

Now, another ten years later, we come back to this subject because we feel that this 
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beautiful idea has not been explored enough. We shall replace the EOS used in the early 
works by another one, based on the model of The usual relativistic mean held models 
do not generate equations of state which lead to solitons because they do not contain higher 
derivatives, required in nonlinear differential equations. Therefore we will introduce a new 
interaction term in the Lagrangian of Q]. 

This work is organized as follows. In the next section we present some simple formulas 
of nonrelativistic hydrodynamics and write them in terms of the baryon and energy density. 
This will help to establish the connection between these equations and the relativistic mean 
held EOS. In the subsequent section we discuss the model of jj], introduce the new interaction 
term and, in the mean held approximation, derive the equation of state. In the following 
section, we combine the hydrodynamic equations with the new EOS and derive the KdV 
equation and its analytical solution. The last section contains some concluding remarks. 


( 1 ) 


II. HYDRODYNAMICS 

The basic equations of nonrelativistic hydrodynamics m are Euler equation 

dv fl\ 

and the continuity equation 

^ + V.(piT)=0 (2) 

In these equations, p is the mass density, which, in the nuclear medium, is related to the 
baryon density through 

p = MpB (3) 

where M is the nucleon mass. The enthalpy per nucleon is given by [l^ : 

dh = Tds + Vdp (4) 

where V = 1/pb is the specihc volume. For a perfect huid {ds = 0) the equation above 
becomes dp = pBdh and consequently: 


Vp = pB^h 

Inserting Q and © in m and in m we obtain 


( 5 ) 


( 6 ) 
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and 


(7) 


dpB 

dt 


+ V ■ {pBv) = 0 


The enthalpy per nncleon may also be written as 

h = E + pb 

where E is the energy per nncleon given by: 


III: 

dE 

dpB 


( 8 ) 




(9) 

which inserted in (jS)) leads to: 

Pb 


de 

(10) 


^ = 7r- 

OpB 

and therefore: 

<1 

<1 

to 

(11) 


With these simple manipnlations we can make explicit the dependence of Enler eqnation 
on e^pb), or in other words, on the eqnation of state of nnclear matter. 


As it will be seen later, it is possible to perform a combination of ® and o in a specific 
way to obtain a nonlinear differential eqnation for the baryon density, with soliton solntions. 
The existence or not of these solitons depends on the fnnction e{pB)- A key featnre of 
(P) and (HU is that they contain the gradient of the derivative of the energy density. If £ 
contains a Laplacian of pb, i.e.. 


£ (X ... + ...V^Pb + ... 


( 12 ) 


then eqnation will have a cnbic derivative with respect to the space coordinate, which 
will give rise to the Korteweg-de Vries eqnation for the density. The most popnlar RMF 
models do not have higher derivative terms and, even if they have at the start, these terms 
are nsnally neglected dnring the calcnlations. In the next section we show how to overcome 
these limitations and obtain a slightly more general expression for e{pB) with Laplacians. 


III. MODIFIED QHD 

We start with the well known variant of QHD- 


I P], employed in 


jCo = - 9vV^) - (M - gs(p)]'il’ + 
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- jiv.-F'"' + \mv^V,V^ _ (13) 

where, as usual, the degrees of freedom are the baryon held 'ip, the neutral scalar meson 
held 0 and the neutral vector meson held V^, with the respective couplings and masses. The 
parameters b and c will be taken from As it can be seen, there are no higher order 
derivatives. This is the reason why the traditional QHD and most of its numerous variants 
do not give an equation of state compatible with soliton propagation, like m- For this 
it is necessary, in hrst place, to modify the Lagrangian, introducing a term with Laplacian 
derivative couplings. The simplest one is given by: 

9v 


Cm = 


my 




(14) 


where the subscrit Af denotes “modihcation” in QHD. First of all, this term is designed to 
be small in comparison with the main baryon - vector meson interaction term 
Because of the derivatives, it is of the order of: 


- rsj - 

2 2 
mb mi 


p 


~ 0.12 (15) 

Oy IlOy 

where the Fermi momentum is kp — 0.28 GeV and my — 0.8 GeV. The form chosen for 

Cm is not dictated by any symmetry argument, has no other deep justihcation and is just 

one possible interaction term among many others It is used here as a prototype. In 
the future, we plan to address this problem with a state of the art Lagrangian. As it can be 
immediately guessed, the integration by parts of Cm leads to a D’Alembertian acting on the 
baryon current. We present m with the derivatives acting on just because it can be 
more easily compared to gy'ip'^fiV^'ip it allows for estimates like m- The total Lagrangian 
is then: 

C = Cq + Cm (16) 

Integrating Cm by parts twice, making the standard mean held approximation 


V, V, >= 6,oVo 


X 0 > = 


(17) 

(18) 


and neglecting the time and spatial derivatives of the helds Vq and (po, we arrive at the 
approximate Lagrangian given by: 

C = - gvloVo) - (M - gsfpo)]^^ - 
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( 19 ) 


bcpo^ 

3 




mv^ 


The equations of motion of the helds are: 


- 9vloVo) - (M - 5'50o)]0 = 0 

(20) 

mv'^Vo = gv'^l^i’ 

(21) 

ms‘^(po = 9s^i^ - b(j)o‘^ - c0o^ 

(22) 


where, for simplicity, we have neglected the terms containing Laplacians. This means that 
the perturbation changes slightly the energy of the system, but does not affect the dynamics 
of the environment, i.e., of the unperturbed ifj, 0 and helds. The energy density is 
calculated using the energy-momentum tensor: 


e —< Too >— ~9oo^' + (doi’) 
which, with the help of (USD, (EOl), (EH) is given by: 


dC' 


(23) 


gv^ 2 , ms^\iM*-M) 
^ = - - oPb + 


2my2' 


n 2 


9S 


+ 


7 


( 2 vr)= 


Mf + -^{M* - M) 

Jo 4:g% 


^9l 


9v‘^ „ v 72 „ 9v‘^ „ f d^PB 


(24) 


. 4PB V Pb - jPb 

my 

where ( 7 ) is the baryon spin-isospin degeneracy factor and M* stands for the nucleon effective 
mass: 

M* = M - 5(500 


The two terms in the last line of (j24j) come from Cj^^. Baryon number propagation in nuclear 
matter has been studied in with the help of the diffusion equation: 

^ = OVVb (25) 

where the diffusion constant D was numerically evaluated as a function of density and 
temperature and found to be T ~ 0.35 fm at densities comparable to the equilibrium 
nuclear density and temperatures of the order of 80 MeV. The exact value of T at T = 0 
is unknown but its dependence both with temperature and density is very smooth and an 
extrapolation to T = 0 points to values even smaller than 0.35 fm. This number is small 
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compared to any nuclear size scale and can be interpreted as indicating that density gradients 
do not disappear very rapidly in under matter. Using (ESD twice in the last term of it 
can be rewritten as: 


gv‘^ (d^pB 
—PB 






9v 


■PbD^[W\W^Pb)] 


(26) 


mv‘*' \ dt'^ j my"*' at\ ' j my"" 

which, in the context of the present calculation, can be neglected because V^(V^pb) << 
[V'^Pb). With this last approximation, the hnal form of the energy density is given by 
without the last term. The baryon density up to the Fermi level {kp) is given by: 

likpf 


= PB = 




As usual, the nucleon effective mass is determined by the minimization of e{M*) with respect 
to M*, which results in the self-consistency relation: 


M* = M - 


9s 7 


d-^k- 


M* 


ms^ (27r)3 Jo ^ ^, 2 ^ 


1/2 


9l \h 


ms^ 


c 


„ iM* - MY + ^{M* - MY 
9s 9 s 


The main result of this section is (El, which, in what follows, will be used in the deriva¬ 
tion of the KdV equation. By keeping (but not we take into account the 

small inhomogeneities, which do not change rapidly in time (and therefore we neglect 
Whereas this approximation is quite often used at the nuclear surface, it is not very common 
for the bulk of nuclear matter. We believe, however, that it is quite reasonable. 


IV. THE KDV SOLITONS 

In this section we follow the treatment developed in and 111 to obtain the Korteweg-de 
Vries equation in one dimension by the combination of and dlD- 

Expression (El is valid for any p^. In order to make contact with the nuclear matter 
phenomenology and use the constraints imposed by the saturation curve, we will now use (j24|l 
to compute the energy per baryon E ©• We then Taylor expand it around the equilibrium 
density po up to second order: 


E(ps) = E{p„) + \(^) (PB-P,f 

2. \UpB / PB=P0 


(27) 
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where the hrst order term vanishes because of the saturation condition 


A(^_M 

OPB \PB 


dE\ 
OpbJ , 


= 0 


Pb=Po pb=P0 

The derivative coefficient in (EZD may be rewritten as ji^ : 


d'^E \ Mcs'^ 


PB=PO 


Po 


where M is the nucleon mass and Cg is the speed of sound in nuclear matter. Calculating 
the enthalpy per nucleon (0 by using (EH) leads to: 


Mr 2 

h = Eipo) + ^{SPB^ + Po^ - ApBPo) 


w 


where 


and 




(28) 


(29) 


Vipo) = 


ms^ 


-M) 


2 


9s 


+ 


7 


f'kp' 


{27ry 


£k{P + 


+-^(M* - Mf + -^(M* - MY 
^9s m 


where the integration limit kp (which appears also in M* = M*{kp)) is hxed at the value 
corresponding to pQ. Inserting (1^ in we arrive at: 


dv , ^ ^ 

- + (v.V)v = 


2po^ 
2 


[6pn(Vps) — 4 po(Vpb)] 


9v 


imv 


r)y(vyB) 


(30) 


We now repeat the steps developed in p and [ll| , write (IHn|l and © in one space dimension 
{x) and introduce dimensionless variables for the density and velocity: 


P = 


Pb 

Po 


and 


V = 


We next dehne the “stretched coordinates” ^ and r as 

^ _ ^ 1/2 ~ ^st) 

R 


in Byiii: 


(31) 


(32) 


(33) 
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and 


^3/2 
T = a ' 


Cst 

~R 


( 34 ) 

where R is & size scale and a is a small (0 < a < 1) expansion parameter chosen to be Q: 

ff = hnml ( 36 ) 


where u is the propagation speed of the pertnrbation. After these tedious but straightforward 
substitutions, we expand EH and (jSl around their equilibrium values: 


p = 1 + crpi + ct^P2 + • • • 


V = avi + a V2 + 


(36) 


(37) 


Equations and o will contain power series in a up to Since the coefficients in these 
series are independent we get a set of equations, which, when combined, lead to the KdV 
equation for pi 


dpi dpi 

17 + 


+ 


5'yVo 


\ d^pi 


= 0 


(38) 


.2Mcs‘^mv^R^J di^ 

which has a well known soliton solution. We may rewrite the above equation back in the 
X — t space obtaining a KdV-like equation for pi with the following solitonic solution: 


Pi{x,t) = — — —sech‘^ 


{u — Cs)CsM 


{x — ut) 


9v V 2po 

where pi = api. The solution above is a bump with width A given by: 

9v [ 


(39) 


A = 


2 po 


Since 


and since 


\ {u — Cs)CsM 

0 < sech^{6) < 1 


(40) 


u 


= a <1 


(41) 


we can see that pi is always smaller than one and is consistent with the expansion (jSni)- Re¬ 
lation m implies that u < Umax = 2cs. This upper bound is arbitrary, being a consequence 
of the choice dsni), but it must exist. As mentioned in the introduction, the perturbation of 
the medium may be caused by the passage of a nucleon projectile and we expect that for 
higher projectile energies the modihcation of medium will be so strong that our perturbative 
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approach will break down. The upper bound in the pulse velocity implies a limit for the 
momentum of the incident nucleon given by: 


Pmax — 


2csM 


(42) 


VI - 4c5^ 

As it can be seen from (P|). the soliton width depends on parameters of the microscopic 
dynamics, as the coupling and mass which come from the Lagrangian (uni). The sound 
velocity comes from the EOS derived from the same Lagrangian. The expression PI|) also 
reveals the supersonic nature of this phenomenon, since it is only possible for u> Cg. 

In the case where the density pulse is caused by a nucleon traversing a nucleus, if the 
width is comparable to the nucleon size and much smaller than the nuclear radius, the pulse 
will be localized and the impinging nucleon may emerge from the nucleus. If, on the other 
hand, A is large, of the order of the nuclear radius, then the incident energy will most likely 
be absorbed by the nucleus as a whole, which will undergo some collective excitation and 
decay. 

In the numerical evaluation of pn|) we need to know the values of the parameters appear¬ 
ing in the Lagrangians (USD and da. For Cq we consider both the original QHD-I, without 
the nonlinear terms, just as a reference, and the most phenomenologically successful model 
of with its different parameter choices labeled by iVLl, A^L3, iVL3 — II and NL — SH. 
The use of the original QHD-I equation of state i n hy dro dynamical calculations leads to a 
too rapidly expanding system, as it was shown in |l7j |. 

All the input numbers as well as the calculated incompressibilities K, equilibrium densities 
Pq, nucleon effective masses M* and sound velocities Cg are shown in Table 1. 
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QHD 

NLl 

NL3 

NL3 - II 

NL-SH 

K(MeV) 

545 

211 

272 

272 

355 

M(MeV) 

939 

938 

939 

939 

939 

ms{MeV) 

500 

492 

508,2 

507,7 

526 

mv{MeV) 

780 

783 

782,5 

781,9 

783 

9s 

8,7 

10,14 

10,22 

10,2 

10,4 

gv 

11,62 

13,28 

12,87 

12,8 

12,9 

M*/M 

0,56 

0,57 

0,6 

0,59 

0,6 

poifm-^) 

0,19 

0,15 

0,15 

0,15 

0,15 


0 

-12,17 

-10,43 

-10,4 

-6,9 

c 

0 

-36,26 

-28,9 

-28,9 

-15,8 

Cg 

0,25 

0,16 

0,18 

0,18 

0,2 


Table 1: input parameters and results for K, M*, po and Cg. 

It is interesting to observe that the maximal momentum Pmax depends on Cg and therefore 
will be different for different models. 

Choosing values for u, from Cg to Umax, and computing both A and p up to Pmax for 
the different models, we obtain the curves shown in Fig. 1. As it can be seen, the width 
decreases with the projectile energy, as anticipated in |l^. We can also see that, the stiffer 
is the equation of state, the wider is the window for soliton formation. At bombarding 
energies of a few hundred MeV the collisions are probably too violent and the departure 
from the initial equilibrium state is too large for the soliton picture to be true. At energies 
smaller than 200 MeV, the width starts to be larger but still signihcantly smaller than the 
nuclear radius. Moreover, the choice of gv in m was just a hrst guess. The derivative 
coupling could be different and much weaker, leading to narrower solitons. As it appears 
now in (gOD A is directly proportional to gy- If the coupling would be half of the value used 
here, even the softest EOS (NLl) would lead to narrow solitons. After having established 
the connection between KdV solitons and relativistic mean held models, we can test other 
theories and other derivative couplings. 

The observation of KdV solitons in nucleon-nucleus collisions was also discussed in [l^ , 
considering the exisiting experimental results at that time. Nowadays much more data are 
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available and it might be interesting to look at them with the soliton pictnre in mind, 
searching for some nnexpected transparency. 


V. SUMMARY AND CONCLUSIONS 


Assnming that RMF models give a satisfactory description of nnclear matter and hnite 
nnclei and, also, that the nnclens can be treated as a perfect flnid, we have investigated the 
possibility of formation and propagation of Korteweg - de Vries solitons in nnclear matter. 
We have npdated the early estimates of the soliton properties, nsing more reliable eqnations 
of state. In doing so we had to introdnce in the Lagrangian a term with derivative conpling 
between the baryon and the vector meson helds. The term nsed here was the simplest one 
and with it we conld already hnd KdV solitons, which have their properties (in particnlar 
the width) controled by the microscopic parameters of the nnderlying theory, snch as masses 
and conplings. Whereas the qnantitative resnlts are not entirely conclusive, they suggest 
that KdV may indeed exist in nuclear matter. It is fair to say that the subject deserves 
further investigation. Finally we would like to mention that, during the last decade, due 


to t 


20 


le CERN-SPS and RHIC programs, hydrodynamics has experienced a great revival 


2l| and supersonic phenomena 2^ were also brought back to explain recent RHIC 


data. Nonlinear effects, as those discussed here, may be relevant also in this context. 
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FIG. 1: Soliton width as a function of the projectile momentum 
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